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Introduction 

Many years ago in paper [I] some special set-theoretical operation Ak 
was introduced. It was used to reduce manipulations with phase prefactors 
in the Fock space ( Grassmann algebra) to pure set theoretical ones. More 
detailed discussion of this operation was given in [2j, Appendix A. The main 
purpose of this paper is to show that the operation mentioned makes it 
easy to obtain very complicated analytic expressions in effective and uniform 
manner. In addition, this technique helps to reveal hidden symmetry of the 
electronic Fock space. Fock space equipped with new algebraic structure can 
be considered as a natural domain for a certain class of Post Hartree-Fock 
(HF) methods. 

In Section I we discuss properties of the operation Ak both for the case 
of molecular spin orbital (MSO) basis and molecular orbital (MO) basis. 

Section II is dedicated to classic configuration interaction (CI) matrix 
elements. It does not contain any new results and is included only as sim- 
ple but important application of our technique of manipulations with phase 
prefactors. 

In Section III excitation operators and coupled cluster (CC) approach jl] 
- PO] are discussed. General analytic expression for CI expansion of CC wave 
function is given. For the CC CI coefficients recurrence relations are derived. 
Analytic expressions for derivatives of CC CI coefficients with respect to CC 
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amplitudes are given. CC density operator of the first order necessary for 
calculation of molecular properties (see, e.g., 0) is constructed. 

Sections IV and V are dedicated to new algebraic structure on the finite- 
dimensional electronic Fock spaces. 

In Section VI and Appendix some hints concerning computer implemen- 
tation of our approach are given. 

Basic Definitions 



Let N = {1, 2, . . . , n} be the molecular spin-orbital (MSO) index set. On 
the set V(N) of all subsets of N let us consider the operation (symmetric 
difference) 

RAS = (RUS)\(Rf]S) (1) 

where R, S G V(N). This operation endows V(N) with Abelian group struc- 
ture with empty set as its unit. Each element of this group is of order 2 
(RAR — 0). 

Let us define the operation Ak as 

A K = A {1,2,..., A;} (2) 

k£K 

where K C N. For example, if K = {2, 4, 5, 7} then 

A K = {1, 2}A{1, 2, 3, 4}A{1, 2, 3, 4, 5}A{1, 2, 3, 4, 5, 6, 7} = {3, 4, 6, 7}. 
The mapping 

tp : K -> A K (3) 

is a group homomorphism. Indeed, A0 = and {A K )A(A L ) = A K & L . 
In general case, for K = k\ < k<i < . . . < k s 



f rjn 

2) 



A K = { 



U {&2i-i + !,-••, &2i} if s is even 

'31 w 

\J{k 2 i + l,...,k 2 i+i} if s is odd 



Directly from the definition of operation (1) the following useful relation 
may be obtained 

\KHR\ + \KnS\ = \Kn(RAS) I (mod 2) (5) 
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For a fixed basis set of n orthonormal spin-orbitals the corresponding 
finite-dimensional Fock space jFjy is spanned by determinants \R) where R 
runs over all subsets of the spin-orbital index set N. Its p— electron sector 
J~N, P is spanned by determinants \R) with \R\ = p. Basis determinants will 
be labelled by subsets and all sign conventions connected with their represen- 
tation as the Grassmann product of ordered spin-orbitals will be included in 
the definition of the creation- annihilation operators. 

Creation-annihilation operators associated with spin-orbital index % are 
defined by the following relations 

a]\R) = (1 - ( i:R )(-iy\RU {i}) (6a) 

a t \R) = ( i>R (-iy\R\{i}) (6b) 

where 

e=|{l,2,...,i-l}ni2| (7) 

is the sign counter and 

«/,«>={; i\i R R <8) 

is the well-known combinatorial ( function of partially ordered by inclu- 
sion set V(N). Note that by an abuse of notation we use symbol ((i,R) 
instead of (({i},R). It is pertinent to mention as well that \R\ = \S\ and 
C(R, S) = 1 imply R = S. 

The main purpose of this section is to get explicit and easily evaluated 
expression for the action of products of the creation-annihilation operators 
on arbitrary determinant wave function. We start with calculation of the 
sum \K fl A.R | + \R fl Aft-|. From Eq.(5) it readily follows that for arbitrary 
K,RcN 

\KnA R \ + \RnA K \= IW 

fl A{ r }| + \ {r} fl A{ fe }| (mod 2) (9) 

(k,r)€KxR 

But 

|{Jb}nA {r} | + |{r}nA w | = {i l\l r r (10) 

As a result, 

\KC] Aft| + \RD A K \ = \K\\R\ + \KDR\ (mod 2) (11) 
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Let us return to the creation-annihilation operators and their products. 
If % £ R then 

\{l,2,...,i- 1}DR\ = \R(~)Ai\ = \R\ + \{i}n A R \ (mod 2) (12a) 
where Eq.(ll) was used. If % e R then 

|{l,2,...,i-l}ni2| = |(i2\{0)nAi| (mod 2) 
With the aid of Eq.(ll) we can rewrite the last expression as 

|{i, 2, . . . ,i - 1} n r\ = \R\{i}\ + \{i} n A RX{t} \ 

= \R\ + \{i} fl Ar\ (mod 2) (126) 
Now formulas (6a) and (6b) may be recast as 

a\\R) = (-1)1*1(1 - C^)(-l) IWnAfll |i?U {i}) (13a) 

a,\R) = (-l)l^ C ^(-l)l« nA «l|i?\{z}) (136) 

Are Eqs.(13a) and (13b) more convenient than Eqs.(6a) and (6b)? In our 
opinion in general situation it is more a matter of taste than convenience. 
But if the reference determinant \R) is fixed and it is necessary to calculate 
the action of different creation- annihilation operators on this determinant, 
then it is possible to find at first the bit vector representing subset A# and 
then each time check only one bit of this vector in contrast to multiple checks 
if Eq.(7) is used. 

We are interested in expressions of the type a\ k . . . a| x a Jfe . . . | R) involv- 
ing particle number preserving products of the creation-annihilation opera- 
tors. Successive application of Eqs.(13b) and (13a) gives 

4 k ■ ■ ■ al a j k ■ ■ ■ a h\ R ) = 

k k 

= II(i-Uw-i)IIU^i x (-m(RVk) u i k ) (14) 

1=1 1=1 

where 
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I 

and 

k 

e = [\R\J1-1\ + \{R\J k ) u + n A MJi _J 

+IO*}nA (iAW _ 1 |] (15) 

Let us try to simplify the expression for e. It is clear that the right-hand 
side of Eq.(14) is nonzero only if | J{\ = I and \Ii\ — I for any I — 1, ... ,k. As 
a result, 

k 

[\R\J1-1\ + \(RVk) U /,_!!] = k (mod 2) 

i=i 

and the expression for e takes the form 

^ i—i k ^ 

^+EEi^ nA >i+EEiw nA >i 

Z=l r=l 1=1 r=l 

k l-l k k 

+ E E i w n + E I0i> n A fl | + e I w n A «i ( 16 ) 

/=1 r=l Z=l Z=l 

Note that the -R-dependent part of this formula is very easy to handle: it 
is sufficient to construct the set Ar and then just to check if the indices 
ji, . . . ,jk and ii, . . . , i k belong to this set. The remainder terms in this for- 
mula depend only on the order relation between indices of the creation- 
annihilation operators involved. 

Expression (16) takes more simple and compact form if ji < . . . < j k and 
%i < . . . < i k . Indeed, in this case, as easily follows from Eq.(4), subsets 
{ji} flAj w and n A^ are empty. As a result, 

e = k + \i k nA Jk \ + \i k nA R \ + \J k nA R \ (17) 

If k = 1 then 

a\a 3 \R) = (1 - CMAO } )0,*(-l)1(^\{j}) U {<}> (18) 

where 

e = l + |{i}nA {j - } | + |{i,j}nA fl | (19) 
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If k = 2 then 

a\a]aia k \R) = 

= (1 - a(fl\{M})u{i})(l - Q,R\{k A Ki,R\{k}(kA-V £ \(R\{k, i}) u {«, j}) (20) 
where 

e= |{i}nA {i} | + |{i,j,z}nA {fc} | + |{i,j}nA {0 | 

+\{i,j,k,i}nA R \ (21) 

In Handy orbital representation [3| each MSO index set is split in its a 
and (3 components. In particular, in expression K n A R subsets K, R go to 
(K a , Kp) and (i? a , Rp), respectively. Let us introduce the sets Kp = m + Kp 
and Rp = m + Rp where m is the number of molecular orbitals (MO). We 
have K = K a AKp and R = R a ARp and 

\K(1 A R \ = \K a n A Ra \ + \K a n A Rg \ + \Kp nA Ra \ + \Kp nA Rg \ (mod 2). 

It is clear that \Kp D A Ra \ = and \K a fl A R | = 1^11^1 (mod 2). As a 
result, 

\KHA R \ = \K a nA Ra \ + \KpnA R!3 \ + \K a \\Rp\{mod2). (22) 

Standard Matrix Elements 

One-electron operator is of the form 

h=^2<i\h\j> a\ aj (23) 

i,j£N 

and its action on determinant \R) may be presented as 

h\R) = -J2 E (-i) IWnA «> l+l{i ' i}nAfl| < > l(iAW) u {0) 

jeii i6(Af\fi)U{j} 

For matrix element < S\h\R > we have the expression 
<S\h\R>=-J2 E (-l) |{i}nA «> l+l ^ }nABl < > C(*W»u Wl s 

j6ii i€(5\fl)u{j} 

(24) 
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Matrix element < i\h\j > is zero if MSO involved are of different spin. As 
a result, in MO basis only cases j G R a ,i G S a and j G Rp,i G Sp are to 
be treated. Note that by an abuse of notation we use the same letters for 
indices of MSOs and MOs. 

Case W . From Eq.(22) it follows that 

\({ki}, 0) n A fl | = \{h} nA Ra \ + \Rp\(mod 2) (25)) 

and 

e aa = |{i}nA - } | + \{i,j}n A Ra \ (26) 
Case '(3j3\ We have 

|(0,{fci})nA fl | = |{A: 1 }nA^|(mod2) (27) 

and 

e m = ({OnA^I + l^^nA^I (28) 

Thus, in MO basis Eq.(24) becomes 

< (S a ,Sp)\h\(R a ,Rp) >= 

jeR a ie(S a \R a )u{j} 

-<R a ,S a J2 E (- 1 ) 6 " 3 < W > C(^\{,})u«,5, (29) 
Two-electron operator is of the form 

g= ^2 < ij\kl > a\a]aia k (30) 

i,j,k,leN 

and its matrix element < S\g\R > may be written as 

< S\g\R >=J2 ( _1 ) £ < yl fcZ > C(iA{M})u{*j},s (31) 

fc.ieft »,je(S\H)u{fc,i} 

where £ is given by Eq.(21). 

Two-electron integrals < ij\kl > equal to zero if MSOs with indices i, k 
or with indices j, I are of different spins. In MO basis two-electron integrals 
are written as (ik\jl) and there are four cases to be analyzed. 
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Case , aaaa\ All orbitals are of the same spin a and in this case 

\({h, k 2 }, 0) n A R \ = \{ki, h} nA R J + 2|^|(mod 2) (32) 

and 

Saaaa = |{0 H A {k} \ + \{j} fl A {fc} | + |{«} fl A {fc} | + \{j} fl A {i} | 

+|{z} n A w | + |{i} n A m | + \{i,j} nA R J + |{fc, z} n A Ra \ (33) 

Case 'a.af3(3\ Orbitals with indices i, k are of a spin and with indices j, I 
are of (3 spin. We have 

|({fci}, {h}) n A R \ = nA R J + |{A; 2 } nA R(J | + I^Kmod 2) (34) 

and 

Saapfs — H A {fc} | + |{j} H A w | + | {i, A;} n A K J + I {j, /} n A R , | (35) 

Case , P(3aa\ Orbitals with indices i, k are of j3 spin and with indices j, I 
are of a spin. We have 

£/3/3aa = | {0 D A {fe} | + |{j} H A {J} | + | {i, fc} n A fl/3 1 + | {j, 1} n A Ra I (36) 

Case All orbitals are of the same spin (5 and in this case 

|(0, {h, k 2 }) n A fl | = \{h, k 2 } n A^|(mod 2) (37) 

and 

e pm = \{l} n A w | + |0} n A w | + \{i} n A w | + \{j} n A w | 
+|{*} n A {0 | + \{i} n A {j} \ + nA fi J + \{k, 1} n A^| (38) 

In MO basis the expression (31) takes the form 

< (S^S^lgK^Rp) >= 

Cr^s, Yl (- 1 ) £aaaa (^b'0C(ii«\{fc,o)u{ij},5 a 

k,lGRa i,je(Sa\Ra)U{k,l} 
(k^l) (ijtj) 

+ J2 Yl (- 1 ) eaaW (* fc iiOC(fla\{fc})u{i} ) s a C(fl /J \{i})u{,}^ 

fce-Ra ie(S a \fl a )U{fc} 

3e(s /3 \fl /3 )u{i} 
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+ J2 E (- 1 ) £,3 " aa (^b'0C(fl«\{o)uo} ! 5 a C(^\{fe})u W ,5, 

kSRp ig(S /3 \fl /3 )U{fe} 

ieR a je(s a \R a )u{i) 

+(R a ,s a E (-1) £ ^%'0C(^\ { M}M^^ (39) 

MSflg i,je(S g \Rp)U{k,l} 

To reduce the dimension of configuration interaction (CI) space one may 
turn to the notion of the active space. Let us suppose that each determinant 
under consideration is presented in the form 

(R a ,R p ) = (lUR a a ,lUR%) (40) 

where / = {1,2,..., n^} is the set of the so-called inactive MO indices that 
occur in each determinant with the occupancy 2. 

From Eq.(29) it easily follows that either i, j e / or i,j £ I. As a result, 

<(lUS a a ,lUS$\h\(lUR a a ,lUR a p ) >= 

Ck^Cr^J^ 2 < ^ > 

iei 

E E ( _1 ) £QQ < W > C(fls\{j})u{<},ss 

-Ci?g,5 2 E E (- 1 ) 6 "" < W > C(^\0})u{<},S| (41) 
ieflg ie(S%\R a g )u{j} 

In the same manner Eq.(39) may be treated. We have 

< (I U S a a , I U S$\g\(I U R a a , I U R%) >= 

ijei 

+Cr^ E E (-l) £aaaa (^b'OC(ii S \{fe,0)u{ij},5 3 

fe.iefig i,je(sg,\R%,)u{k,l} 

+ J2 E (- 1 ) eaaW ( ifc b'OC(fla\{fc})u{i},s s C(^\{i})u{j}.^ 

fceftg ie(Sg\fl°)u{fc} 
ie-R^ je(S<»\.R<*)u{i} 
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+ E E (- 1 ) e " /,aa ( ifc b'OC(fla\{i})uo-},ssC(flg\{fc})u{i},s| 

ieflS ie(sg\flg)u{!} 

+Cfla,sa E E (-i) £ ^%'OC ( ^\ { m})u { ^,^ (42) 



Excitation Operators 

The operator performing /c-fold excitation from some fixed single deter- 
minant reference state \R > is defined as 

a J k = a l---4 1 a j k --- a n (43) 
where = { j± < . . . < jt} C R and I k — {h < ■ ■ ■ < ik} C N\R. We have 
a \\R >= (_l)*+l^nA Jfc |+|/ fc nA fl |+|j t nA il ||( i2 y fc ) y Jfe > 

Product 

gives non- vanishing result when applied to |i? > if and only if families 
I kl , . . . , 7 fc and J kl , . . . , J^ M include mutually disjoint subsets. After simple 
manipulations we get 

a^: I Z\R>=(-ir\(R\J)Ul> (45) 



where 

fi fi 

j=l i=l i>j 

(46) 

and I = (J I ki , J= U Jfc.. 

i=l i=l 

The expression for e can be simplified if we suppose that the reference de- 
terminant \R) is built on the first p spin-orbitals that is \R) = |{1, 2, . . . ,p}). 
We, however, do not presuppose any special choice of indices of occupied HF 
MSOs in our further discussion. 
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Definition 1. Family { Jj}f =1 is a set partition of J into fx blocks if 

(1) Ji ^ for all z = 1,2,...,//; 

(2) ^ D Jj = for all % ^ j; 

(3) J = U J<. 

i=l 

Normally set partition has its blocks listed in increasing order of small- 
est element in each block. We will call such partitions block ordered ones, 
reserving the term 'set partitions' for arbitrary families with properties (1) - 
(3) of Definition 1. 

If | J\ = k then the number of all possible block ordered set partitions of 
J is equal to the kth Bell number Bk (see, e.g., [TT]). 

Definition 2. Family {( Jj, Jj)}f =1 ) is a consistent set partition of a pair 
of subsets (J, /) if 

(1) U}?= 

x is a block ordered set partition of J; 

(2) {/«}f =1 is a set partition of /; 

(3) \ Ji\ = for alii = 1, 2, . . . , //. 

Excitation operators appear explicitly in coupled cluster theory jlj - [TU] 
where the exponential anzatz of wave function is used: 

\qf>=e x \R>. (47) 

Here 

i (*) 

fc=l J fc CH 
I k CN\R 

is the excitation operator associated with chosen reference determinant 
\R>, and I < p is the maximal excitation order. 
Let us consider the standard expansion 

eX = E^ ( 49 ) 

and try to present the expression for the restriction of X M on one-dimensional 
reference subspace C\R > in a form convenient for further analysis. We have 

min(p,fil) (k) (fci + ...k ll =k) 

*"\r>= E E E *t-^ a t^i i2> (50) 

k=n JCR J k ,...,J k CJ 
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where the inner sum goes over families that determine consistent set parti- 
tions of subsets J and I with 1 < ki < I. Summation of both sides of the 
last equality over /i gives 



p min(p,fil) (k) (&i+...fc M =fc) 

i«>+e E E^ E ^-^"tij 



ICN \ R ^,-,^c/ 

(51) 

The right-hand side of Eq.(51) contains similar terms. Indeed, if \i is fixed 
and some family {{Jki}i=n {-^fci}f=i) is selected then any permutation cr <E 
generates new family ({</fc CT(i) }f=i, {h a(i) }i=i)- Of ^' families thus obtained 
let us choose one and put 

^ = Ei( / ^) UJ ^)) nA (^, u ^,) 1 (52) 

i>j 

Then the summation over /i! families on the right-hand side of Eq.(51) is 
reduced to calculation of the following sum 



s,{{J ki }U {I ki }U = £ (53) 

where 



fJ " aes„ 



£ = ^2ki + ^2\i k nA Jk \ + |/n A R \ + |Jn A R \ (54) 
i=i i=i 

But from Eq.(ll) it easily follows that 

\(I ki U J ki ) n A (4 . UJfc .)| + |(/ fc . u j fc .) n A (/fciUJfci) | = 0(mod 2). 

As a result, J fci }f =1 , {/fci}f=i) = ( — 1) £ and final expression for e x |i? > is 

v (*) 

e x | J R>=^ £ T l k (J,I)\(R\J)Ul>, (55) 



where 



fc=0 Jcfi 

ICN\R 



(ki+...kfj,=k) 
'fci 

12 



and 

7^(0,0) = 1 (566) 
In Eq.(56) the index set M k is denned as 

Mk = {/i : [/i, min(p, fil)] 3 k} 

It is easy to see that k is the maximal element of Its minimal element 
may be found from the inequality 



k = l- 



+ r < fj, • I 



where < r < I, and is given by the following relation 

_f[f] if k = 0(mod I) 
^™ n ~\[f] + l itk^0(modl) 

This fimin corresponds to the minimal possible number of blocks of sizes 
not greater than / (see Appendix). Thus, the summation in Eq.(56a) goes 
actually over all consistent set partitions of the pair (J, I) such that each 
block size does not exceed I, and e is defined by Eq.(46). 

At this stage it seems reasonable to introduce new CC amplitudes that 
differ from the initial ones defined by Eq.(48), only by phase prefactors: 

t i _> (_-Q|J|+|/nAj|+|/nA fl |+|JnA,R| t / ^ 

By an abuse of notation for phase modified CC amplitudes we will use the 
same symbols as for the initial ones. With new definition of CC amplitudes 
Eqs.(55) and (56) stay unchanged except for e in Eq.(56) that becomes 

£ = ^|(/,,u4)nA (W | (59) 

i>j 

The expression (55) establishes a connection between CC and CI expan- 
sions in most general case. It is pertinent to mention that in CC expansion all 
determinants are involved, even if the CC operator (48) is of low excitation 
order. 

Now we are ready to get the expression for the first order density operator 
p C R C = ^[e x \R><R\e x ], (60) 
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where the contraction operator is defined by 

^~W< 5 I = - \ -(-l) 5 \j)(i\ if S = (R\{j}) U {i} ( 61 ) 

if \Sr\R\ <p- 1 

with 

<J = |{i}nA {j - } | + |{i,j}nA fl | 

and i ^ j. 

After simple manipulations we come to the following expression: 

P (fe) 

Mf = E E Kw j )] 2 E b->c/i 

fc=o jcr jeR\Jui 

ICN\R J S 

P (*) 

E E T ^ J ' J ) E(- 1 ) 71+72 ^(( J \{0)uo-},/)ii)(ii 

fc=l ./C-R je-R\J 
ICN\R ie.J 

p-1 (k) 

-EE T ^ J ) E (-i) 7i+72 ^ + i(^u{j},/u«)|j)(.| 

fc=0 Jcfl jefl\J 

ICN\R ieN\R\I 
V (*) 

E E ^(^ / )E(- 1 ) 71+72T ^i( J \{o,AO})b)^i 

fe=l JC-R JS-T 
ICN\R ieJ 

V (fc) 

"EE T ^ J ) E (-^^(^(AOWU^IjXil (62a) 

fc=l Jcfl je/ 

ICN\R ieN\R\I 

where 

7l = |{*} n A m |, 72 = |{i, j} n A {RV)UI \ (626) 
or, after collecting similar terms, 

jeR fe=o jcr\{j} 

ICN\R 
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i,jeR 



p (k) (k— 1) 

+ E woiEE E KM/uo-})] 2 

j'GA r \i? fc=l JCRICN\R\{j} 

' V (fc-i) (*) 

E E E (-i) 7i+72 ^u{o,/)^(ju{j},/) 

fe=l JcR\{m} ICN\R 



- E i^'X*i 



iGN\R 



p-1 (fc) 

E E (-ir^(j,/)^ +1 (ju{j},/u«) 



fc=0 JCH\{j} 
/CiV\H\{i} 



+ E 



iSJV\H 



P 

E E (-ir^(Ju{<},/uo})^-iU/) 



fc=l jck\{0 

JCJV\H\{j} 



+ E i^"><*i 



l,j£N\R 



p (k) (fc-1) 

EE E (-i) 7i+72 ^(^ / u {j})t^( j, / u «) 

fc=l JCRlcN\R\{i,j} 



(63) 

In the analogous manner the expression for the second order CC density 
operator may be obtained. 

For further analysis it is convenient to introduce the coefficients 



[nMn= E 



(fci+...+fc M =fc) 



(/* + ?)! 



E (-iK:;-^ 



(64) 



where // m j n is defined by Eq.(58) and r is the vector (matrix) of CC ampli- 
tudes. It is clear that [T^ (r)] J7 = T^(J, J), and that [7o ? (t)] J7 = ^. 
Let us present the coefficient [^^(t)] 3I in the form 



min(k,l) (fc„) ^ 

= E E 
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I*CN\R 



d(Xtj, 
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(65) 



Calculation of the derivatives of [^^(t)] n with respect to (phase mod- 
ified) CC amplitudes we start with the analysis of phase prefactor. To this 
end let us recast Eq.(59), using Eq.(ll), as 

e= J2 l(4 ; u J fci )nA (4 . UJ( 



--kj)\ 



i>3 



+ 1(4, u J ki ) n A^yjajl + 1(4,, u J k J n A ( ^ UJO |. 

i 

Taking into account that 



f/fc. uiJnA 



2fc. + 1 



= k it (mod 2) 



we can write the following elegant expression 



d 



■j- \tLJXt)] TT d\ 
d ( Xt jJ 



9+1 V J)j\J,I\h 



if J D J* and J D I* 
if J 7$ J* or J 7$ i* 



where 



r)(J*,h) = fc* + | (J* U 2*) n A (Ju j)|. 
Now Eq.(65) may be recast in the form of the recurrence relation: 



(66) 
(67) 



J7 



min(k,l) (fc») 



fc» = l i.CJ 

/*c/ 



^+i(^W*a^ (68) 



CC CI coefficients correspond to the case q = in this equation. The recur- 
sion can go either to the vacuum (k — k* = 0) or to small values of k — k* for 
which the required CC CI coefficients can be evaluated directly. 
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Structure of Commutative Algebra on 
p - Electron Sector of the Fock Space 



Let us denote by the symbol A P R the p— electron sector Tn, p of the Fock 
space where some determinant \R) is selected. In more formal terms A P R is a 
pointed space that is the pair (J~n, p , \R))- The selected point \R) is referred 
to as either the HF vacuum or the HF reference state. 

A P R is spanned by the basis vectors 

e^R) = \(R\J) U /) (69) 

labelled by subsets J C R,IC N\R with |J| = |/|. Note that ej(R) is 
nothing more than a special convenient notation for the basis determinant in 

If S = (R\J) U / then it is easy to ascertain that 

4(S) = efX^%(R) (70) 

where 

J[ = J' n (R\J), 4 = J'nl (7ia) 

and 

i[ = i'n J, I' 2 = I'D (N\R\I). (71b) 

In particular, e^(S) = ej(R) and ej(S) = e®(R). 

Therefore, arbitrary finite linear combinations of vectors corresponding 
to different vacuum states are allowed since with the aid of Eqs.(70) - (71) 
such combinations are easily reduced to some common vacuum state. 

We furthermore assume that the vacuum vector \R) is fixed and, by an 
abuse of notation, for basis vectors (69) we will use the symbol ej. 

A P R is Hermitean space with scalar product 

Vj (72a) 



i 

2 

(72b) 
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v (fc) 



and norm 



\x\\ = 



k=o JCR 

ICN\R 



p (fc) 

E E W 

k=0 JCR 

ICN\R 



Subspaces A l R of the dimension 

«=e(;)(v> 

k=0 v 7 v 7 

spanned by the basis vectors with 1 < | J\ — \I\ < Z, form the following 
chain 

A° R C ^ C . . . C ^ (74a) 

where 

^ = CeJ. (746) 

For each pair (k, I) with k < I we introduce subspace Wjj spanned by 
vectors ej with < | J\ = \I\ < I. It is clear that W R is the orthogonal 
complement to A k R l in ^4 Z R and 

p 

A P R = W [ R Kk) (75) 

fc=0 

For a fixed CC excitation level I the subspace W R is 'the subspace of CC 
amplitudes 

Let us define the multiplication in A P R by putting 

r v _{ (-i)l(^wl e ^; if j n j> = and / n I' = f76) 
3 J '~ \o if Jnj'^0 or /ni^r 1 j 

It is easy to see that with such a multiplication the vector space A P R becomes 
a commutative and associative algebra with ejjj as its identity (compare with 
algebras studied by Paldus 6J). Commutativity follows directly from Eq.(ll), 
associativity from the relation 

\{J X U h) n A (j2U/2 )| + |(Ji U J 2 U /a U / a ) n A (J3U/3) | 

= 1(^2 u h) n A (J 3 U/3 )| + |(Ji u A) n A ( j 2U j 3 u7 2 u73)l(^ orf 2). 

T/ie structure introduced depends on the choice of the HF vacuum vec- 
tor and for a fixed MSO basis set there are exactly (™) such (isomorphic) 
structures. 

As has already been mentioned, for vector space structure on A P R the 
change of vacuum vector simply means change of notation for the same basis 
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determinants (see Eqs.(70) - (71)). In contrast, the algebra structure defined 
by Eq.(76), is deeply connected with the basic vacuum vector. For example, 
e®(S) = \S) is the idempotent of algebra A p s whereas the same determinant 
\S) = ej(R) in algebra A P R is its nilpotent basis element. Isomorphism be- 
tween A p s and A P R is induced, for example, by a permutation a of the MSO 
index set N such that cr(S) = R. 
The relation 

where k\ + k 2 < p, shows that A P R is a graded algebra. The subspace W R is 
its maximal nilpotent ideal and the algebra under discussion is just a direct 
sum of the field of complex numbers and this nilpotent ideal 

A P R = Ce§ © (78) 

Note also that A P R is algebra with involution induced by the complex conju- 
gation. 

It is well-known [T2j that in algebras of such type it is easy to define 
algebraically both the exponential mapping and its inverse. Namely, 

p T k 

k=0 

_^ J , v\;(l,i>) 



is the exponential mapping 

-> eg + Wr ■ and 

log-.el + r^^-^T (80) 
k=i 

is the logarithmic mapping ejjj + W R ' P ^ — > W R . Here 

T k = T*. . .*T . (81) 

k 

For mappings (79) and (80) all classic relations hold true: 

exp(ri) * exp(r 2 ) = exp(ri + r 2 ), (82a) 

[exp(T)Y l = exp(-r), (82b) 
exp(log(e® + T)) = T, (82c) 
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and 

log(exp(r)) = r, (82d) 

where r, Ti, r 2 are arbitrary elements from W R ' P \ 

If r G W£' is a vector of CC amplitudes then it is easy to show that the 
CI coefficients (56) are just the components of the vector exp(r) in algebra 

v (k) 

exp(r) = J2 E T i( J ^)^ ( 83 «) 



A p ■ 



k=0 JCR 

ICN\R 

where 

T l k {J,I) = {exp{r)\e I J ). (836) 

In fact, exp{r) in the last equation for each fixed k may be replaced by the 
sum of the first k terms in its expansion. Eqs.(82) - (83) may be considered 
as a source of different relations involving CC CI coefficients. 

Thus, the Fock space equipped with the non-trivial structure of commuta- 
tive algebra, defined by Eq.(76), can be considered as a natural domain for 
the CC and related approaches. 

Spectral theory in algebra A P R is very simple. Indeed, let us write arbitrary 
element x G A P R as 

0i 

X — -|- / 

where r G Y\? R l,p \ Its spectrum is defined as the set of all A G C such that 
(Aeg — x) is not invertible in A P R . It is easy to see that 

(Ae» -*)-' = £ — T * =R(*,A) (84) 

exists for all A ^ x\. Thus, spectrum of element x G A P R is 

Spx = {(x\el)} (85) 

The mapping C\{x0} — > A P R defined by Eq.(84) is called the resolvent of 
element x G A P R . 

From Eq.(85) it follows that if element x G A R corresponds to normalized 
p— electron wave functions then its spectrum belongs to the unit disk of the 
complex plane : 

Sp x C D = {z G C : \z\ < 1} 
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In general case the structure of normed space on A P R is not compatible 
with the algebra structure defined by Eq.(76), that is the inequality 

||x ★ y|| < ||x|| • ||y|| (86) 

not necessarily holds true. For example, if n = 2 and p = 1, then for x = 
eg + ej^j the inequality (86) is violated. Indeed, in this case x*x — + 2e|^| 

and = a/5 whereas ||x|| 2 = 2. Our hypothesis is 

Hypothesis. For sufficiently large n inequality (86) is satisfied for any 
x, y G A P R . 

If this hypothesis is true then for values of n of actual interest A P R is 
a commutative unital Banach algebra(see, e.g., [HJ QSj). In this case dif- 
ferential calculus of full value could be developed on A P R . It would be too 
restrictive to consider only finite expansions of the type of Eqs.(79) - (80) 
and it would be possible to define the exponential and logarithmic mappings 
by their series and prove classic relations (82a) - (82b) for arbitrary x G A P R . 
And even more, Banach algebras are very closely related to algebras of ana- 
lytic functions of complex variables and this also could be of interest 
for applications. 

The basic idea of CC and related methods is to parametrize subsets of A P R 
using elements('amplitudes') from W R where I is the maximal excitation 
level. General parametrization of such a type may be defined as an injective 
mapping 

vr : - A P R (87) 

such that element 7r(r) is invertible for any r G YV R . The last requirement 
guarantees that (7r(r)|ejjj) 0, or, in other words, that the HF reference state 
appears with nonzero coefficient in expansion of 7r(r) for each r. It is clear 
that the exponential mapping (79) is a parametrization (see Eq.(82b)). It 
possesses the following characteristic property: for any excitation level I the 
set exp(W R ' ) is a multiplicative Abelian group. This property is referred to 
as 'size consistency' and is usually considered as essential advantage of the 
exponential parametrization over other ones. 

Another parametrization used in the quadratic configuration interaction 
method (QCI) [13^ is given by 

^(r, 1) = f> fc = (e§ - r)- 1 (88a) 

k=0 
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where 

Q l k (J,I) = (n(T,l)\e I J ). (886) 
are the corresponding CI coefficients. The mapping inverse to 1Z(t, 1) is 

n- l :el + T^Y,(-l) k - l r k (89) 
k=i 

The exponential and QCI parametizations are just special cases of general 
polynomial parametrization 

v 

P a .. T ^J2a k r k , (90) 

fc=0 

where ) G C p+1 . Element P a {r) is invertible for any r and 

the mapping (90) is injective if and only if do • a.\ ^ 0. The inverse to P a may 
be written in the form 

P- 1 :a el + r^Y.{-l) k - 1 b k r\ (91) 
k=i 

where coefficients b k are easily obtained from the relation 

iti-^hW-aoeff^T. (92) 
k=i 

We have h = ±,b 2 = ^, etc. 

It is pertinent to mention that if the excitation level I is less than the 
number of electrons p then any parametrization covers only a part of the 
full CI space. If x is some CI vector and (7r _1 (a:)|e5) = for all J, I with 
\J\ = \I\> I then 7T — 1 (a?) G VV^ 1,0 , that is this CI vector x is covered by the 
parametrization used. 

The main attractive feature of any polynomial parametrization consists 
in possibility to get approximation to the lowest eigenvector of Hamiltonian 
by iterative solution of non-linear system of relatively small order 

(e I J \H\P a (r)) = E(e I J \P a (r)), (93) 

where 1 < I J| = 1/1 < /. 
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In algebra A P R it is easy to derive relations involving vector r and energy 
value E and different from the standard CC type equations. For example, 
let us suppose that QCI parametrization (88) covers some eigenvector of the 
Hamiltonian H, that is there exists r e wj ? 1,/ ' ) such that 



HH(t, 1) = EK(r, 1) 
It is easy to ascertain that the mapping 



(94) 



k=l 



k 



(95) 



is inverse to the mapping r — > ETZ(t, 1). Since, by assumption, Eq.(94) holds 
true, we can apply the transformation (95) to the both sides of Eq.(94) to 



k=i 



T. 



(96) 



This relation may serve as a base for Newton-Raphson type iterative algo- 
rithms for solution of Eq.(94). 

In algebra A P R the expression for the first order density operator for arbi- 
trary polynomial parametrization takes the form: 



PP C R ° = 

jeR 



v (*) 

E E 14X41 

k=0 JCR\{]} 
ICN\R 



\P a (r))\j)(j\ 



+ E 

j£N\R 



V (*) (fc-i) 

EE E 

k=l JcRIcN\R\{j} 



'J l\ e J 



\P.(r))\m 



T 



i,j£R 
(iltj) 



(fc-1) (fc) 

E E E (-i)" 2 i^u W )(eSu { , } i 

fc=l JCR\{i,j} ICN\R 



\Pa(T))\j){i\ 



i€N\R 



p-1 (fc) 

E E (-i) T ' rt i4x4 u u gi 



JCJV\H\{»} 



l^(r))|j)(i| 
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+ E ( p ^\ 



j€N\R 



fc=l JCfl\M 
JCJV\H\{j} 



l^(r))|j)(z| 



+ E 



i,jeN\R 



P (*) (*-!) 

EE E (-i) 72 i4 uw )(4 uw i 

fc=l JcRlcN\R\{i,j} 



\P a (r))\j)(i\, 

(97) 

where the expressions for 71 and 72 are given by Eq.(62a). 

If / and g are differentiable mappings M. — > ^ then it is easy to see that 



A[/(A) * ff (A)] = [^/(A)] *a(A) + /(A) * [^(A)]. 



(98) 



In particular, for rea/ vector r of CC amplitudes the following relations hold 
true 

_r* = At*" 1 * eS, (99a) 



dA 



and 



T exp(r) = exp(r) * ej, 
j 

exp(Ar) = exp(Ar) *r, 



(996) 
(99c) 

(99d) 



In concluding this section let us note that it is easy to define algebra 
^\rr) corresponding to orbital representation of p— electron sector of the 
Fock space spanned by determinants with fixed spin projection Ms = \{p a ~ 
pp). The multiplication in A^^^ is 

(i a ,i p ) {i' a ,i'p) _ 

(-!) CT e (j„uj>.j»uj>,) if ^HJ CT = and 7^4 = 



if J CT nj;^0 or J CT n/;^ 



(100) 



where a — a,f3. 
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Pointed Fock Spaces as Associative 
Superalgebras 



Let us describe possible extension of the multiplication rule (76) to the 
full Fock space. To this end we consider a pointed Fock space Ar = (J~n, \R)) 
with the basis ej(R) = (R\J)Ul) where J C R and / C N\R( on the whole 
2 P ■ 2 n ~ p = 2™ basis vectors) and present Ar as a direct sum 

A R = A + R @A R (101) 

where subspaces .4.^ and A R are spanned by basis vectors ej(R) with even 
and odd values of | J U I\, respectively. Using the same definition (76) for 
the product of arbitrary basis vectors, we come to associative algebra with 
identity. From Eq.(ll) if follows that Ar is a skew- commutative (or super- 
commutative) algebra: 

e^R) *e5(i2) = (-if^'^e^R) * e^R). (102) 

We have 

^ "^Ri *^R ^ *^R "^^R^ *^R ^ "^R "^^R^ *^~R *^R %/ ^'R^ (103^ 

which means that Ar is Z 2 — graded (super) algebra and that .4^ is a subal- 
gebra of Ar. 

Let us consider the pointed Fock space A® that is the Fock space where 
the 'absolute' vacuum state \R) = |0) is selected. In this case we have 

4m^) = {^' nA ^W ;^n/; = (104) 

It is interesting to compare the star product in A% with the standard 
exterior (wedge) product in the Grassmann algebra Tn- We have 

|/)A|rH {(-«ur> -jnr-j (105) 

where e(I, I') is the number of pairs (i, i') e / x /' such that i > i! . It is easy 
to show (see Eq.(ll)) that 

£(/,/') = |/'n Aj| = |/n A r \ + |/| • \I'\(mod 2) 
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and this means that even in the case of algebra A% the star product is differ- 
ent from the standard exterior product in the Grassmann algebra T^. For 
example, if / = {1, 2, 7} and V = {4, 5, 6} then A// = {1,2,3,4, 6}. We have 

\inA r \ = 2 

5(/,r) = |{(7,4),(7,5),(7,6)}| = 3. 

Another possible extension of the multiplication rule (76) is given by the 
relation 

j r _{ (-l)l^'^') nA (-)le^; if J n J' = and I n I' = (m) 

Both definition (76) and definition (106) lead to the same structure on 
the p— electron sector of the Fock space. Their extension to the full Fock 
space gives, however, different results. In particular, with the definition 
(106) algebra A$ is identical to the Grassmann algebra. 



Computer Implementation of Star Product 

For two arbitrary vectors x, y G A P R we have 

P (k) 

x * y = x\y + y\x + 

k=l JCR 

ICN\R 

(107) 

Methods of calculation of this rather complicated product are based on 
two fundamental combinatorial notions: listing and ranking. Listing means 
generation of combinatorial objects (subsets of a given index set in our case) 
in some fixed ordering. Ranking algorithms allow to set up arrays indexed 
by subsets (vectors of algebra A P R in our case). In other words, rank of a 
subset is the position this subset occupies in a given ordering. 

There exists simple formula for calculation of subset rank under lexical 
ordering of the set of all k— element subsets of the index set iV (see, e.g., 
HH1): 

te * W= 6)-§^*-|/n"4>l + i)- (108) 



g^(_ 1)fcl+ l(^u/0nA (JUI) | x / > A/ ;i 

fc 1= l JiCJ 
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In algebra A P R components of arbitrary vector x are indexed by pairs of 
subsets (7, 7) with 7 C R and 7 C N\R. Let us put 

lex R (J) = lex T - HR) (r-\j)) (109a) 

and 

lex N \ R (I) = /ex s -i (A r\ jR) (s _1 (7)) (1096) 

where R — r\ < . . . < r p and N\R — s\ < . . . < s n - p . Vector x G A P R has its 
components indexed as 

(7,7) - (lex R (J) - l)(^~ P ^+lex N \ R {I). (110) 

In our opinion ranking and listing algorithms required for evaluation of prod- 
ucts x-ky should not necessarily be closely related. For example, we can use 
Gray codes generation algorithms [To*! ITT] and lexical indices (110). 

Conclusion 



Advanced technique of manipulations with phase prefactors arising in 
many problems of quantum chemistry is demonstrated. On the base of this 
technique three methods of evaluation of CC CI coefficients and CC density 
matrices are suggested. 

The first method (see Eq.(56)) requires rather complicated combinato- 
rial algorithms for set partitions generation and permutations of multisets 
(see Appendix) . The success of its implementation strongly depends on the 
program code efficiency. 

The second method (see Eq.(68)) is based on the recurrence relations for 
CC CI coefficients and its efficiency is closely related both to the efficiency of 
recursive algorithms in programming language used and available hardware. 

The third method employs new algebraic structure revealed on the Fock 
space. The efficiency of implementation of a rather special product of two 
vectors determines the overall efficiency of any algorithms using this multi- 
plication rule. 

The first two methods do not require full size CC CI vector to be kept 
in fast memory but they may be rather time consuming. The third method 
can be programmed in a very efficient way but at least two CC CI vectors 
should be stored. 
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Appendix: RG strings and set partitions 



All information contained in this Appendix is well-known and is given 
here only to show how direct evaluation of expressions of the type of Eq.(56) 
can be performed in non-recursive way. Standard references are \W\ I17j. A 
lot of interesting information may also be found on Frank Ruskey site |18j . 

Let { Ji}i =1 be a set partition of the set iV = {1,2, ... ,n} having its blocks 
listed in accordance to the following rule: block Jj contains the minimal 
element of the set A^\(|J*^ 1 1 Jj). In particular, J\ 3 1. In this Appendix 
the term 'set partition' means 'block ordered set partition' unless otherwise 
stated. 

Definition. Restricted Growth (RG) string of length n is an integral 
vector a = (ai, 02, • • • , a n ) such that 

ai < dj < max {a,-} + 1 (1) 

l<j<i 

where a\ — l. 

With each RG string it is convenient to associate the following vector 

v = (yi, . . . ,v n ) (2) 

where Vi = max{ai, . . . , a^}. 

Lemma. There exists a bijection between the set of all RG strings of 
length n and the set of all set partitions of N. 

Proof. See P5||IZ|. 

Non-recursive generation of all RG strings in lexical order can be described 
as follows. We start with the minimal RG string 

a min = (1, 1, . . . , 1) (3) 

Let us suppose that on some step we have current RG string 

(4) 

where \J am \ > 1 and |Ja m+ J = ... = |J 0n | = 1. It is easy to see that in 
RG string (4) a m+1 > v m . Indeed, in the opposite case the smallest element 
of J am+1 would be not greater than m but this is in contradiction with the 
assumption that J am+1 is a block of size 1 containing index m+ 1. As a result, 
the current RG string (4) takes actually the form 

(l,...,a m ,u m + l,...,u m + n-m) (5) 
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The corresponding set partition with v m + n — m blocks is 

{ Ji, . . . , J am , ■ ■ ■ , Jv m , {m + 1}, . . . , {n}} 

The next (in lexical order) RG string is 

(l,...,a m + l,l,...,l) (6) 

In going from RG string (5) to RG string (6) the size of the block J\ increases 
by n — m and the current number of blocks decreases by the same value and 
becomes equal to v m if a m ^ v m . If a m = v m then RG string (6) corresponds 
to the set partition with i/ m + 1 blocks. 

Generation of RG strings should be terminated when the number of cur- 
rent blocks becomes equal to n that corresponds to the maximal RG string 

O-max = (1,2, ... ,n) (7) 

and the set partition with exactly n blocks Jj = {i}. 

We are interested in generation of set partitions with the additional re- 
striction on block sizes. Namely, for any partition of iV the size of each its 
block should not be greater than the maximal excitation order I. This im- 
plies, in particular, that with such a restriction the number of blocks can not 
be less than 

f [s] if n = (mod I) 

\ [f ] + 1 if n £ 0(mod I) l8j 
With lexical ordering of admissible RG strings the minimal one is 

(l,. l ..,l,2,. l ..,2,...,k,. l .^k) iin^0(modl) ( 9 ) 

l l n-l-[f] 

and the maximal one is given by Eq.(7). 

Let us suppose that the current RG string is of the form of Eq.(5) and 
that the corresponding vector of block sizes is 

S=(|Jl|,...,|Ja m |,...,|^ m |,l,...,l). (10) 

Since each block size should not exceed /, it is necessary to scan vector (10) 
from the position a m + 1 till v m to find the smallest index a m such that 
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Ja m \ < /. If such an index is not found then we put a m = v m + 1. The next 



where K\ = min{l — \ Ji\,n — m},/t 2 = min{l — | «^| , — m — Ki}, etc. If 
after sorting blocks J 1} . . . , J Um {Jv m +i) there are still k = n — m — ^Kj>0 



unfilled positions in the string under consideration, then we should fill them 
by Si = min{l, k} integers v m + \{y m + 2), etc. 

The next algorithm required is the algorithm of generation of all set par- 
titions (not necessarily block ordered ones) corresponding to a given vector 
s — (si, . . . , S/j) of block sizes. But this problem is equivalent to the problem 
of generation of multiset permutations and the corresponding non-recursive 
lexical algorithm is easily obtained by modification of, say, algorithm 5.1 
from the book |19j . 
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